@-PLOS | ove 


CrossMark 


click for updates 


o OPEN ACCESS 


Citation: Yang G, Xu Z, Lu W, Li X, Sun C, Guo J, et 
al. (2015) Quantitative Analysis of Differential 
Proteome Expression in Bladder Cancer vs. Normal 
Bladder Cells Using SILAC Method. PLoS ONE 
10(7): e0134727. doi:10.1371/journal.pone.0134727 


Editor: Jon M. Jacobs, Pacific Northwest National 
Laboratory, UNITED STATES 


Received: January 20, 2015 
Accepted: July 13, 2015 
Published: July 31, 2015 


Copyright: © 2015 Yang et al. This is an open 
access article distributed under the terms of the 
Creative Commons Attribution License, which permits 
unrestricted use, distribution, and reproduction in any 
medium, provided the original author and source are 
credited. 


Data Availability Statement: All relevant data are 
within the paper files and its Supporting Information 
files. 


Funding: This work was supported by the National 
Science Foundation for Young Scientists of China 
(No. 81402115, 81201572), the Natural Science 
Foundation of Jiangsu Province, China (No. 
BK20140172) and the 111 Project of China (No.111- 
2-06). The funders had no role in study design, data 
collection and analysis, decision to publish, or 
preparation of the manuscript. 


PLOS ONE | DOI:10.1371/journal.pone.0134727 July 31, 2015 


Quantitative Analysis of Differential 
Proteome Expression in Bladder Cancer vs. 
Normal Bladder Cells Using SILAC Method 


Ganglong Yang’, Zhipeng Xu”, Wei Lu’, Xiang Li?, Chengwen Sun“, Jia Guo’, Peng Xue°*, 
Feng Guan'* 


1 The Key Laboratory of Carbohydrate Chemistry & Biotechnology, Ministry of Education, School of 
Biotechnology, Jiangnan University, Wuxi, China, 2 Shaanxi Provincial People’s Hospital, Xian, China, 

3 Wuxi Medical School, Jiangnan University, Wuxi, China, 4 Department of Urology, Affiliated Hospital of 
Jiangnan University, Wuxi, China, 5 Laboratory of Proteomics, Institute of Biophysics, Chinese Academy of 
Sciences, Beijing, China 


* xuepeng@moon.ibp.ac.cn (PX); fengguan @jiangnan.edu.cn (FG) 


Abstract 


The best way to increase patient survival rate is to identify patients who are likely to progress 
to muscle-invasive or metastatic disease upfront and treat them more aggressively. The 
human cell lines HCV29 (normal bladder epithelia), KK47 (low grade nonmuscle invasive 
bladder cancer, NMIBC), and YTS1 (metastatic bladder cancer) have been widely used in 
studies of molecular mechanisms and cell signaling during bladder cancer (BC) progression. 
However, little attention has been paid to global quantitative proteome analysis of these 
three cell lines. We labeled HCV29, KK47, and YTS1 cells by the SILAC method using three 
stable isotopes each of arginine and lysine. Labeled proteins were analyzed by 2D ultrahigh- 
resolution liquid chromatography LTQ Orbitrap mass spectrometry. Among 3721 unique 
identified and annotated proteins in KK47 and YTS1 cells, 36 were significantly upregulated 
and 74 were significantly downregulated with >95% confidence. Differential expression of 
these proteins was confirmed by western blotting, quantitative RT-PCR, and cell staining 
with specific antibodies. Gene ontology (GO) term and pathway analysis indicated that the 
differentially regulated proteins were involved in DNA replication and molecular transport, 
cell growth and proliferation, cellular movement, immune cell trafficking, and cell death and 
survival. These proteins and the advanced proteome techniques described here will be use- 
ful for further elucidation of molecular mechanisms in BC and other types of cancer. 


Introduction 


Bladder cancer (BC) is the fifth most common type of human cancer. There were an estimated 
74,690 newly diagnosed cases and 15,580 deaths from this disease in the United States in 2013 
[1]. Of total BC patients, >70% have nonmuscle-invasive disease and ~25% present initially 
with muscle invasion. Patients with the muscle-invasive form have a 50% risk of distant 
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that no competing interests exist. reason for the worldwide prevalence of BC [3]. The majority (90%) of BCs are classified histo- 
logically as urothelial carcinomas (UCs), derived from the bladder urothelium [4]. Bladder epi- 
thelial tissues have a clear hierarchical organization consisting of three morphologically 
distinct cell types: basal, intermediate, and umbrella cells, corresponding respectively to early, 
middle, and late differentiation states [5]. Malignant transformation may occur in each of these 
cell types, resulting in a diversity of tumor phenotypes [6]. According to the latest report of the 
American Cancer Society, the relative 5-year survival rate for BC with early detection (stage I, 
(T1, NO, MO)) is ~88% [7]. Therefore, identification of novel early-stage molecular markers is 
desirable for improved risk stratification. 

Candidate biomarkers for BC detection evaluated to date include telomerase, bladder tumor 
antigen (BTA), nuclear matrix protein 22 (NMP-22), and fibrin degradation product (FDP). 
The reliability of tests based on these biomarkers is poor because of low sensitivity and high 
false-positive rates [8-11]. Proteins can potentially be identified specific to aggressive or nonag- 
gressive types of cancer. Proteome analysis is challenging because of the limited amount of 
available clinical sample [12]. Monitoring of the proteome of BC cells could provide additional 
information for clinical diagnostic purposes. 

Recent advances in mass spectrometry (MS) for protein identification and quantification 
facilitate in-depth analysis of large numbers of proteins, and have been used for examination of 
the whole proteome in several systems. Such methods include 2D difference gel electrophoresis 
(2D DIGE), the similar iTRAQ (isobaric tag for relative and absolute quantitation), isotope- 
coded affinity tagging (ICAT), and stable isotope labeling by amino acids in cell culture (SILAC) 
[13-15]. In comparison with peptide-based absolute quantitation methods, SILAC has the 
advantages of mixing samples at the very beginning, and reduced sample-to-sample variability. 
Metabolic labeling with stable isotopes has been described as the "gold standard" in protein 
quantification [16]. Arginine (Arg) and lysine (Lys) are the stable isotope-labeled amino acids 
most frequently used in SILAC-based studies, because subsequent trypsin digestion of isolated 
proteins (which cleaves at basic residues) for MS analysis generates peptides with a single labeled 
amino acid, simplifying analysis and quantification [17]. In the present study, three stable iso- 
topes each of Arg (RO, R6, R10) and Lys (KO, K4, K8)in three separate cultures (“light” (L), 
“medium” (M), and “heavy” (H)) were used to analyze proteome differences at various stages of 
BC. Distinctive L, M, and H forms of each peptide as detected by MS reflected relative amounts 
of the corresponding protein in three isotopically encoded BC cell stages. 

Three human cell lines were studied: normal bladder epithelial HCV29, low grade nonmus- 
cle invasive bladder cancer (NMIBC) KK47, and metastatic muscle invasive bladder cancer 
YTS1. Each of the three cell lines was cultured in media added with three combinations of unla- 
beled Lys and Arg ("light"), D4-Lys and 1 C,-Arg ("medium"), and 2 C6""N2-Lys and BC Ne 
Arg ("heavy"). Proteins with >98% label incorporation were analyzed and quantified by 
2D-HPLC-LTQ Orbitrap MS (Fig 1). Differential expression of the identified proteins, which 
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Fig 1. Schematic procedure for quantitative analysis of proteins in BC cells vs. normal bladder cells. 


doi:10.1371/journal.pone.0134727.g001 
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are presumably related to BC development, was confirmed by western blotting, quantitative 
RT-PCR, and cell staining with specific antibodies. 


Material and Methods 
Cell culture 


HCV29, KK47, and YTS1 cells were established as described previously [18-20] and kindly 
donated by Dr. Sen-itiroh Hakomori (The Biomembrane Institute; Seattle, WA, USA). Cells 
were cultured in RPMI 1640 medium supplemented with 10% FBS and 1% penicillin/ strepto- 
mycin at 37°C in 5% CO, atmosphere. For SILAC labeling, cells were cultured in SILAC- 
labeled RPMI 1640 with 10% dialyzed FBS and 1% penicillin/ streptomycin containing 

“light” (KORO), “medium”(K4R6), or “heavy” (K8R10) Lys and Arg. To prevent Arg-to-Pro con- 
version, L-Pro (200 mg/L) was added to the medium as described previously [21]. Cells were 
cultured for at least 5 passages to eliminate nonlabeled Lys and Arg. 


Cell lysis and protein extraction 


Total proteins of the three cell lines were lysed and extracted using T-PER Reagent (Thermo 
Scientific; San Jose, CA, USA) according to the manufacturer's instruction. In brief, cells 
(~1x10’) were detached with trypsin, washed twice with ice-cold 1xPBS (0.01 M phosphate 
buffer containing 0.15 M NaCl, pH 7.4), lysed with 1 mL T-PER reagent containing protease 
inhibitors (1 mM PMSF and 0.1% aprotinin), incubated for 30 min on ice, homogenized, and 
centrifuged at 12,000 rpm for 15 min. The supernatant was harvested and stored at -80°C. Pro- 
tein concentration was determined by BCA assay (Beyotime; Haimen, China). 


SDS-PAGE and in-gel digestion 


Proteins were separated by 10% SDS-PAGE, visualized by Coomassie staining for 2 hr, and 
destained overnight. Excised gel slices were washed with 25 mM ammonium bicarbonate/ 50% 
acetonitrile (ACN). Dried pieces were incubated with 20 uL of 10 mM dithiothreitol (DTT) at 
56°C for 1 hr and then with an equal volume of 20 mM iodoacetamide (IAM) at room temper- 
ature in the dark. Gel slices were washed and trypsinized with 20 uL (10 ng/uL) trypsin (Pro- 
mega; Madison, WI, USA) for 30 min at 4°C. Excess trypsin solution was removed, 20 uL of 25 
mM NH,HCOs was added, and the mixture was incubated overnight at 37°C. Extracted prod- 
ucts were dried with a SpeedVac concentrator (CentriVap Cold Trap, Labconco; Kansas, MO, 
USA) [22]. 


MALDI-TOF/TOF-MS 


Dried samples were dissolved with 0.1% trifluoroacetic acid (TFA), spotted onto an MTP 
AnchorChip sample target, and air-dried. Peptides were recrystallized with matrix o-cyanohy- 
droxycinnamic acid (CHCA; 1 uL of 10 mg/mL) and characterized by MALDI-TOF/TOF-MS 
(UltrafleXtreme, Bruker Daltonics; Bremen, Germany). Ionization was achieved by irradiation 
with a nitrogen laser (A = 337 nm) operating at 20 Hz. Mass spectra were acquired using the 
FlexControl and FlexAnalysis software programs. 


In-solution digestion 


Proteins from three types of stable isotope-labeled cells were mixed at 1:1:1, reduced, and alkyl- 
ated by incubation with equal amounts of 10 mM DTT and 20 mM IAM. Alkylated proteins 
were digested by trypsin added at a ratio of 1:50 (w/w) and incubated overnight at 37°C [23]. 
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Total peptides were concentrated and desalted using a 10KD size-exclusion spin ultrafiltration 
unit and dried using a SpeedVac concentrator. 


LC-MS/MS analysis 


2D-LC-MS was performed using LTQ Orbitrap MS (Thermo Fisher Scientific; Waltham, MA, 
USA) as described previously [24]. Digested peptides (100 ug) were injected into a biphasic 
capillary column (i.d. 200 um) packed with Cg resin (ReproSil-Pur, 5 um, Dr. Maisch GmbH) 
and strong cation-exchange resin (Luna 5 um SCX 100A, Phenomenex). Peptide effluents from 
the biphasic column at each step were directed into a 15-cm Cg analytical column (i.d. 75 um, 
ReproSil-Pur, 3 um) at flow rate 500 nl/min. Nano-ESI was performed with spray voltage 2.0 
kV and heated capillary temperature 200°C. One full MS scan (300-1800) in the Orbitrap was 
followed by five MS/MS scans of the five most intense ions selected from the MS spectrum in 
LTQ. Charge state screening was enabled for +2, +3, +4, and above [25]. 


Data analysis 


Raw MS data were analyzed using the MaxQuant software program (V. 1.2.2.5) [26,27]. A false 
discovery rate (FDR) of 0.01 for proteins and peptides and a minimum peptide length of 6 
amino acids were required. MS/MS spectra were searched by Andromeda [28] against the IPI 
human database (V. 3.85). The MaxQuant program determined the SILAC state of peptides 
from mass differences between SILAC peptide pairs, and this information was used to perform 
searches with fixed Arg6 and Lys4 or Arg10 and Lys8 modifications, as appropriate. Quantifi- 
cation in MaxQuant was performed as described previously [26]. 

Differential regulation within each experimental M/L ("medium/ light") ratio and H/L 
("heavy/ light") ratio of the identified proteins was normalized using z-score analysis, as 
described previously [29,30]. In brief, M/L and H/L ratios were converted into log2 space, and 
average ratios and SD (standard deviations) were calculated for each data set. The log2 M/L 
and H/L ratio of each protein were converted into a z-score, using the following formula: 


z—scores(a) of [b] 


log, *“ " )[b] — Average of (log,of each number, a...n) 


Standard deviation of (log,of each number, a...n) 


where b were deemed as a single protein in a data set population (a. . ..n). The z-score was a 
measure of how many SD units (0) of the log2 M/L or H/L ratio of the protein was away from 
the population mean. A z-score > 1.9600 represented that differential expression of the protein 
lied outside the 95% confidence interval, a score >2.5760 represented expression outside the 
99% confidence interval, and a score >3.291o represented 99.9% confidence. Z-scores >1.9600 
were considered to be significant [29]. 


Functional annotation and Ingenuity Pathways Analysis 


Identified proteins were further analyzed using the SWISS-PROT database to classify their bio- 
logical process, cellular component, and molecular function [31]. Significant over-represented 
gene ontology (GO) terms were identified using the Database for Annotation, Visualization and 
Integrated Discovery (DAVID) gene bioinformatic resources [32, 33]. Proteins determined to be 
differentially regulated as described in the preceding section were tabulated in Excel and their 
International Protein Index (IPI) numbers were uploaded into DAVID (http://david.abcc.ncifcrf. 
gov/home.jsp) for functional annotation analysis. Data sets containing gene identifiers and corre- 
sponding expression values were then uploaded into the Ingenuity Systems application. Each IPI 
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number was mapped to its corresponding gene object in the Ingenuity Pathways Knowledge 
Base. Networks of the proteins were generated algorithmically based on their connectivity. Fish- 
er's exact test was used to calculate a p-value indicating the probability that a particular biological 
function and/or disease assigned to that network was due to chance alone. 


Western blotting 


Western blotting was performed as described previously [34]. In brief, proteins were separated 
on 10% SDS-PAGE gel, transferred onto a PVDF membrane, and the membrane was blocked 
using 5% nonfat milk in TBST and probed using specific antibodies. The primary antibodies 
were rabbit anti-insulin-like growth factor 2 mRNA-binding protein 1 (IGF2BP1) 
(#AP10466b, Abgent; San Diego, CA, USA), rabbit anti-melanoma-associated antigen 4 
(MAGEA4) (#AP6166a, Abgent), rabbit anti-Thy-1 membrane glycoprotein (THY1) 
(#AP2050a, Abgent), 14-3-3 protein sigma (SFN) (#sc-365539, Santa Cruz Biotechnology, Dal- 
las, TX, USA), rabbit anti-fibronectin (FN1) (#F3648, Sigma-Aldrich), mouse anti-vimentin 
(VIM) (#V5255, Sigma-Aldrich), CD70 antigen (CD70) (#sc-7681, Santa Cruz Biotechnology), 
and mouse anti-B-catenin (CTNNB1) (#610153, BD Biosciences; San Jose, CA, USA). The sec- 
ondary antibodies were appropriate horseradish peroxidase (HRP)-conjugated rabbit anti- 
mouse or goat anti-rabbit (#40216 and #A0208, Beyotime). Bands were visualized using an 
enhanced chemiluminescence detection kit (Westar Nova, Cyanagen; Bologna, Italy). 


Quantitative real-time PCR 


Cells (1x10° per well in a 6-well plate) were cultured and treated as described above. Total 
RNA was isolated using an RNApure Tissue Kit (CWBiotech; Beijing, China) according to the 
manufacturer’s instructions. Primers were designed using the DNAMAN program (V. 6.0.3; 
Lynnon Biosoft, Canada). First-strand cDNA was synthesized from total RNA using ReverTra 
Ace-a (Toyobo; Osaka, Japan). Quantitative real-time PCR was performed by LightCycler- 
based SYBR Green I dye detection with UltraSYBR Mixture (CWBiotech). Gene expression 
was quantified by the 2 44CT method [35]. 


Cell staining 


Cells were cultured on sterilized coverslips in 24-well plates until 70-80% confluence, washed, 
immobilized, permeabilized with 0.2% Triton X-100 for 10 min at room temperature, and 
blocked with 5% nonfat milk overnight at 4°C. Fixed cells were incubated with diluted primary 
antibodies for 12 hr, incubated with FITC-labeled secondary antibodies at 4°C for 6 hr in the 
dark, washed, stained with 4 ug/mL DAPI at room temp for 10 min, washed with 1xPBS, and 
photographed with a fluorescence microscope (Eclipse E600, Nikon; Tokyo, Japan). 


Results 
Determination of isotope incorporation efficiency 


To analyze dynamic changes in BC oncogenesis at the proteome level, the SILAC method was 
applied to three bladder cell lines to obtain labeled cell populations. Sufficient labeling is a pre- 
requisite for reliable quantification using this method. In the case of incomplete labeling of pro- 
teins, particularly for labeling efficiency <95%, quantitation of low-abundance proteins would 
be masked by contaminated “light” peptides. To determine incorporation efficiency of labeled 
Lys and Arg, “light” (HCV29), “medium” (KK47), and “heavy” (YTS1) proteins were sepa- 
rated, and the high-abundance protein was in-gel digested. MALDI-TOF/TOF-MS results for 
peptide GVVDSEDLPINISR in heat shock protein 90 (P08238) indicated that complete 
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Fig 2. Mass spectrometric analysis of stable isotope-labeled proteins (SILAC method). (A) 
Determination of incorporation efficiency by MALDI-TOF/TOF-MS. Peaks annotated as RO (left), R6 (middle), 
and R10 (right) are peptide GVVDSEDLPLNISR of heat shock protein 90 from HCV29, KK47, and YTS1 
cells. (B) Identification and quantification of proteome in BC cells by 2D-HPLC LTQ Orbitrap MS. Peaks 
annotated as KO, K4, and K8 (left) and RO, R6, and R10 (right) are doubly charged peptide 
VNQIGSVTESLQACK of alpha enolase and GGPEVQQVPAGER of fatty acid synthase. 


doi:10.1371/journal.pone.0134727.g002 


incorporation of isotopically labeled Arg and Lys was achieved in KK47 and YTS1, and no 
Arg-to-Pro conversion occurred (Fig 2A). LC-ESI-MS/MS analysis of the doubly charged pep- 
tide VNQIGSVTESLQACK of alpha-enolase (P06733) and GGPEVQQVPAGER of fatty acid 
synthase (P49327) showed that these doublets of actual peak clusters were from HCV29, 
KK47, and YTS] cells, respectively (Fig 2B). 


SILAC cell model for quantification of proteome in BC progression 


Proteins isolated from the three cell lines were mixed (1:1:1) and digested using a 10 KD filter 
(Millipore; Billerica, MA, USA). Peptides were analyzed by ultrahigh-resolution liquid chroma- 
tography-tandem MS (nLC-ESI-MS/MS) on a hybrid linear ion trap LTQ Orbitrap instrument. 
A total of 3721 unique proteins were identified in two independent replicate experiments (Fig 
3A and S1 Table). Of these, 1766 proteins (47.5% of the total) that were identified in both 
experiments and satisfied the criteria established for protein quantitation were subjected to fur- 
ther bioinformatic analysis. The distribution histograms of log ratios for both M/L and H/L fit 
a Gaussian distribution. Most of the identified proteins were within the +1 range of log ratios 
(Fig 3C and 3D). Using 1 as the threshold log ratio, expression of 255 proteins was higher in 
both KK47 and HCV29, and expression of 434 proteins was lower in both KK47 and HCV29 
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Fig 3. Distributions of proteins identified in various experiments described in the text. (A) Venn diagrams of numbers of identified proteins from 
individual experiments. (B) Ratios of KK47/HCV29 (M/L) and YTS1/HCV29 (H/L) for the set of 1766 proteins. log2 of the SILAC ratio for each protein (n = 2) 
reflects differences in relative expression among KK47, YTS1, and HCV29 cells. (C) Distribution of SILAC M/L ratios. (D) Distribution of SILAC H/L ratios. (E) 
Cluster graph ("heat map") generated by hierarchical clustering of significant regulated proteins after averaging z-scores using a 95% cutoff. 


doi:10.1371/journal.pone.0134727.g003 


(Fig 3B). Population distribution-based z-scores allowed direct comparison of proteins from 
different experiments. Differing confidence level cutoffs were applied to the data by z-score 
analysis to determine which proteins were significantly differentially regulated. The cutoffs 
applied were 95%, 99% and 99.9%, corresponding to z-scores of +1.960, +2.576, and +3.291, 
respectively. Using a 95% cutoff, significant differential regulation was observed for 110 pro- 
teins in KK47 vs. HCV29 (36 upregulated, 74 downregulated) and for 87 proteins in YTS1 vs. 
HCV29 (17 up, 70 down). Differential regulation was observed for 35 proteins in the two BC 
lines vs. HCV29 using a 99% cutoff (2 up, 33 down), but for only five of these proteins using a 
99.9% cutoff (2 up, 3 down) (Table 1). Tables 2 and 3 list the upregulated and downregulated 
proteins determined in the two experiments, with their average SILAC ratios and z-scores. All 
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Table 1. Protein number, log2 ratio mean+SD, and z-scores of SILAC-labeled proteins. 


Cell line Mean log 2 SD Log 2 z-scoresî 

+1.9600 +2.5760 +3.2910 
KK47 (M/L) -0.072 1.237 36, 74 13, 22 7,3 
YTS1 (H/L) -0.151 1.143 17,70 6, 27 1,3 
Both cell lines 2, 33 2,3 0,0 


aThe first and second value shown are, respectively, the number of upregulated and downregulated proteins outside the indicated confidence level. 


doi:10.1371/journal.pone.0134727.t001 


proteins differentially regulated with >95% confidence had a >5-fold alteration of SILAC 
ratios, and most proteins differentially regulated with >99% confidence had a >10-fold alter- 
ation of SILAC ratios. 

Hierarchical clustering analysis of samples was performed to examine correlations of prote- 
ome patterns among the three cell lines. The Cluster graph ("heat map") shows that samples of 
the same cell line cluster together (Fig 3E). Some proteomes were distinctive among the three 
cell lines based on significant alterations in metastatic vs. low grade nonmuscle invasive, 
whereas other proteomes were moderate and consistent. Proteins in the former group may be 
involved in BC development. 


Functional classification and pathway analysis of identified proteins 


Functional interpretation is a crucial step in data analysis when extensive functional annotation 
of the data sets is not available. Taking into account their nonexclusive localization in GO, the 
identified proteins were linked to at least one annotation term each within the GO molecular 
function, biological process, and molecular component categories. The most common molecu- 
lar functions were binding (47.2%), and catalytic activity (30.9%) (Fig 4A). The major biologi- 
cal process categories were cellular (16.3%), single-organism (14.2%), and metabolic (13.8%) 
(Fig 4B). The major cellular component categories were cell (17.6%), cell part (17.6%), and 
organelle (15.8%) (Fig 4C). 

To identify enrichment terms associated with the upregulated and downregulated groups of 
proteins after averaging of z-scores using the 95% cutoff, lists of proteins were uploaded to the 
DAVID website using the complete human proteome as background. To help clarify which 
molecular functions and biological processes were most affected during BC maturation, over- 
represented GO terms were identified based on threshold count > 2 and Expression Analysis 
Systematic Explorer (EASE) < 0.1. The over-represented molecular functions, biological pro- 
cesses, and cellular components in the significant enriched GO terms of upregulated proteins 
were analyzed. The most highly ranked molecular function was neutral amino acid transmem- 
brane transporter activity (2 proteins). The most highly ranked biological processes were cellular 
amino acid derivative metabolic process (4 proteins), ribonucleoprotein complex biogenesis (4 
proteins), response to extracellular stimulus (4 proteins), and nitrogen compound biosynthetic 
process (4 proteins). The most highly ranked cellular components were intracellular non-mem- 
brane-bounded organelle (14 proteins) and non-membrane-bounded organelle (14 proteins). 

Next, over-represented molecular functions, biological processes, and cellular components 
in the significant enriched GO terms of downregulated proteins were analyzed. The most 
highly ranked molecular functions were calcium ion binding (16 proteins), structural molecule 
activity (11 proteins), and identical protein binding (10 proteins). The most highly ranked bio- 
logical processes were cell adhesion (14 proteins), biological adhesion (14 proteins), response 
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Table 2. Upregulated proteins in BC cells with >95% confidence’. 


Swiss-prot Gene name Protein name M/L H/L log2 M/L log2 H/L z-scores z-scores 
average average average average M/L H/L 
C9JGI3_HUMAN TYMP Thymidine phosphorylase 39.51 1.43 5.30 0.51 4.35 0.47 
FOLR1_HUMAN FOLR1 Folate receptor alpha 31.73 2.91 4.99 1.54 4.09 1.30 
J3QRJ3_HUMAN THY1 Thy-1 membrane glycoprotein 29.72 2.57 4.89 1.36 4.01 1.16 
RCN3_HUMAN RCN3 Reticulocalbin-3 23.97 3.29 4.58 1.72 3.76 1.45 
STEA4_ HUMAN STEAP4 Metalloreductase STEAP4 23.43 4.23 4.55 2.08 3.74 1.74 
EF1A2_HUMAN EEF1A2 Elongation factor 1-alpha 2 17.58 0.87 4.14 -0.20 3.40 -0.10 
1433S_HUMAN SFN 14-3-3 protein sigma 16.87 3.27 4.08 1.71 3.35 1.44 
DDX21_HUMAN DDxX21 Nucleolar RNA helicase 2 14.33 10.65 3.84 3.41 3.16 2.82 
KCRB_HUMAN CKB Creatine kinase B-type 14.15 1.30 3.82 0.38 3.15 0.37 
CPSM_HUMAN CPS1 Carbamoyl-phosphate synthase 12.22 0.99 3.61 -0.01 2.98 0.05 
[ammonia], mitochondrial 
MTAP_HUMAN MTAP S-methyl-5-thioadenosine phosphorylase 11.66 10.09 3.54 3.33 2.92 2.75 
ASSY_HUMAN ASS1 Argininosuccinate synthase 11.08 1.86 3.47 0.89 2.86 0.78 
HOYDA6_HUMAN NAPRT1 Nicotinate phosphoribosyltransferase 10.23 2.82 3.35 1.49 2.77 1.27 
1A68_HUMAN HLA-A HLA class | histocompatibility antigen, A- 8.34 0.44 3.06 -1.19 2.53 -0.91 
32 alpha chain 
EPIPL_HUMAN EPPK1 Epiplakin 8.12 0.86 3.02 -0.21 2.50 -0.11 
PGH2_HUMAN PTGS2 Prostaglandin G/H synthase 2 m95 1.13 2.99 0.17 2.48 0.20 
HPDL_HUMAN HPDL 4-hydroxyphenylpyruvate dioxygenase- 7.93 1.13 2.99 0.18 2.47 0.20 
like protein 
K1C18_HUMAN KRT18 Keratin, type | cytoskeletal 18 7.51 2.18 2.91 1.13 2.41 0.97 
WDR3_HUMAN WDR3 WD repeat-containing protein 3 7.26 3.36 2.86 1.75 2.37 1.47 
PADI2_ HUMAN PADI2 Protein-Arg deiminase type-2 6.99 1.38 2.81 0.46 2.33 0.43 
EXOS4_HUMAN EXOSC4 Exosome complex component RRP41 6.57 3.20 2.72 1.68 2.26 1.42 
A6PVX1_HUMAN SELENBP1 Selenium-binding protein 1 6.41 1.04 2.68 0.06 2.23 0.10 
KYNU_HUMAN KYNU Kynureninase 6.32 0.35 2.66 -1.50 2.21 -1.16 
4F2_HUMAN SLC3A2 4F2 cell-surface antigen heavy chain 6.24 0.88 2.64 -0.18 2.19 -0.09 
PO210_HUMAN NUP210 Nuclear pore membrane glycoprotein 210 6.16 1.05 2.62 0.07 2.18 0.12 
G3V588_ HUMAN ITPK1 Inositol-tetrakisphosphate 1-kinase 6.09 1.71 2.61 0.77 2.17 0.68 
ANM7_HUMAN PRMT7 Protein Arg N-methyltransferase 7 5.99 2.18 2.58 1.13 2.15 0.97 
S100P_HUMAN S100P Protein S100-P 5.92 0.62 2.57 -0.69 2.13 -0.50 
BOS8I7_HUMAN LAGE3 L antigen family member 3 5.87 4.50 2.55 2.17 2.12 1.81 
J3QKT4_HUMAN PYCRI1 Pyrroline-5-carboxylate reductase; 5.72 2.35 2.52 1.23 2.09 1.05 
MCMBP_HUMAN MCMBP Mini-chromosome maintenance complex- 5.65 4.77 2.50 2.25 2.08 1.88 
binding protein 
PRI2_HUMAN PRIM2 DNA primase large subunit 5.47 4.59 2.45 2.20 2.04 1.84 
K1C17_HUMAN KRT17 Keratin, type | cytoskeletal 17 5.39 0.89 2.43 -0.16 2.02 -0.07 
LAT1_HUMAN SLC7A5 Large neutral amino acids transporter 5.35 1.27 2.42 0.34 2.01 0.34 
small subunit 1 
MYPN_HUMAN MYPN Myopalladin 5.34 0.88 2.42 -0.19 2.01 -0.09 
NPM3_HUMAN NPM3 Nucleoplasmin-3 5.29 3.29 2.40 1.72 2.00 1.45 
UHRF1_HUMAN  UHRF1 E3 ubiquitin-protein ligase UHRF1 4.98 5.13 2.32 2.36 1.93 1.97 
ITA6_HUMAN ITGA6 Integrin alpha-6 4.58 12.65 2.19 3.66 1.83 3.02 
CND3_HUMAN NCAPG Condensin complex subunit 3 4.04 5.67 2.01 2.50 1.69 2.08 
UBS3B_HUMAN UBASH3B _ Ubiquitin-associated and SH3 domain- 3.88 6.30 1.96 2.66 1.64 2.21 
containing protein B 
E9PD53_HUMAN SMC4 Structural maintenance of chromosomes 3.51 5.72 1.81 2.52 1.52 2.09 
protein 
(Continued) 
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Table 2. (Continued) 


Swiss-prot 


KIF4A_HUMAN 
MRP_HUMAN 
CTRO_HUMAN 
CIP2A_HUMAN 
MAGA4_HUMAN 
PLEK2_ HUMAN 
AL1A3_HUMAN 


RGS10_HUMAN 
E41L3_HUMAN 
UCHL1_HUMAN 


Quantitative Proteome Analysis of Bladder Cells 


Gene name Protein name M/L H/L log2 M/L log2 H/L z-scores 
average average average average M/L 
KIF4A  Chromosome-associated kinesin KIF4A 3.26 6.56 1.70 2.71 1.44 
MARCKSL1 MARCKS-related protein 2.74 5.39 1.45 2.43 1.23 
CIT Citron Rho-interacting kinase 2.66 5.31 1.41 2.41 1.20 
KIAA1524 Protein CIP2A 1.87 5.15 0.90 2.36 0.79 
MAGEA4 Melanoma-associated antigen 4 1.48 17.47 0.56 4.13 0.51 
PLEK2 Pleckstrin-2 0.98 5.95 -0.02 2.57 0.04 
ALDH1A3 Aldehyde dehydrogenase family 1 0.93 13.75 -0.11 3.78 -0.03 
member A3 
RGS10 Regulator of G-protein signaling 10 0.79 6.50 -0.34 2.70 -0.22 
EPB41L3 Band 4.1-like protein 3 0.72 6.54 -0.47 2.71 -0.33 
UCHL1 Ubiquitin carboxyl-terminal hydrolase 0.42 9.45 -1.25 3.24 -0.95 


isozyme L1 


*Proteins shown have at least one z-score value (M/L or H/L) >1.9600 in two biological replicates. 


doi:10.1371/journal.pone.0134727.t002 


z-scores 


to wounding (13 proteins), and immune response (13 proteins). The most highly ranked cellu- 
lar components were plasma membrane (40 proteins), plasma membrane part (30 proteins), 
and non-membrane-bounded organelle (24 proteins) (S2 Table). 
Proteins were further analyzed, and metabolic and canonical pathways and interconnect- 
ing proteins were generated, using Ingenuity Pathways Analysis (IngenuityH Systems, www. 
ingenuity.com). The top network functions identified as upregulated proteins in BC cells were 
involved in DNA replication, amino acid metabolism, molecular transport (52 proteins; Fig 
5A), gene expression and hereditary disorders (33 proteins), cell growth and proliferation (20 
proteins; Fig 5B), and post-translational modification and cancer (4 proteins). The top net- 
work functions identified as downregulated proteins in BC cells were cellular movement and 
immune cell trafficking (97 proteins; Fig 5C), lipid metabolism (31 proteins; Fig 5D), cellular 
development, growth and proliferation, and cell death and survival (6 proteins). These find- 
ings indicate that BC cell proteomes were continuously shifting depending on the stage of cell 


metastasis. 


Confirmation of MS results by western blotting 


Variations of the differential proteins described above were confirmed by western blotting. 
THY1, MAGEA4, IGF2BP1 and SFN were detected at higher levels in BC KK47 and YTS1 cells 
than in normal bladder epithelial HCV29 cells, whereas VIM, CTNNB1, FN1 and CD70 were 
detected at lower levels in KK47 and YTS1 than in HCV29 (Fig 6B and 6C). In general, the 
western blotting results were consistent with the variables from MS analysis (Fig 6A; S1 Table 


and S1 Fig). 


Confirmation of SILAC results by qRT-PCR 


The expression of six responding genes at the transcriptional level was evaluated by qRT-PCR. 
In BC KK47 and YTS] cells, expression of MAGEA4, THY1, and IGF2BP1 was significantly 
increased, whereas that of VIM, CTNNB1, and FN1 was greatly reduced (Fig 6D). These find- 
ings are consistent with SILAC results. 
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Table 3. Downregulated proteins in BC cells with >95% confidence’. 


Swiss-prot Gene Protein name M/L H/L log2 M/L log2 H/L z-scores z-scores 
name average average average average M/L H/L 
1C07_HUMAN HLA-C HLA class | histocompatibility antigen 0.17 0.02 -2.53 -6.05 -1.99 -4.83 
A8MUB1_HUMAN TUBA4A Tubulin alpha-4A chain 0.35 0.09 -1.51 -3.42 -1.16 -2.70 
AL1B1_HUMAN  ALDH1B1 Aldehyde dehydrogenase X, mitochondrial 0.11 0.10 -3.14 -3.34 -2.48 -2.64 
ALDR_HUMAN AKR1B1 Aldose reductase 0.11 0.04 -3.13 -4.58 -2.48 -3.65 
ANO10_HUMAN ANO10 Anoctamin-10 0.16 0.17 -2.65 -2.53 -2.09 -1.99 
ARMC9_HUMAN ARMC9 LisH domain-containing protein ARMC9 0.16 0.12 -2.67 -3.10 -2.10 -2.45 
ASC_HUMAN PYCARD _ Apoptosis-associated speck-like protein 0.15 0.38 -2.72 -1.39 -2.14 -1.06 
containing a CARD 
AT2B4_ HUMAN ATP2B4 Plasma membrane calcium-transporting 0.16 0.06 -2.65 -3.96 -2.08 -3.15 
ATPase 4 
B2MG_HUMAN B2M Beta-2-microglobulin 0.42 0.10 -1.24 -3.29 -0.95 -2.60 
BCAT1_HUMAN BCAT1 Branched-chain-amino-acid 0.09 0.57 -3.41 -0.80 -2.70 -0.59 
aminotransferase, cytosolic 
BIN1_HUMAN BIN1 Myc box-dependent-interacting protein 1 0.21 0.12 -2.22 -3.09 -1.74 -2.44 
CATB_HUMAN CTSB Cathepsin B 0.07 0.37 -3.84 -1.45 -3.05 -1.11 
CAV1_HUMAN CAV1 Caveolin-1;Caveolin 0.16 0.38 -2.68 -1.40 -2.11 -1.08 
CBPA4_ HUMAN CPA4 Carboxypeptidase A4 0.11 0.08 -3.18 -3.73 -2.51 -2.96 
CD70_HUMAN CD70 CD70 antigen 0.16 0.10 -2.65 -3.26 -2.09 -2.58 
CD97_HUMAN CD97 CD97 antigen 1.98 0.17 0.99 -2.59 0.86 -2.04 
CD99_HUMAN CD99 CD99 antigen 0.14 0.82 -2.84 -0.28 -2.24 -0.17 
CKAP4_ HUMAN CKAP4 Cytoskeleton-associated protein 4 0.16 0.23 -2.64 -2.12 -2.07 -1.66 
CNN3_HUMAN CNN3 Calponin-3 0.11 0.42 -3.25 -1.25 -2.57 -0.95 
CO3_HUMAN C3 Complement C3 0.84 0.15 -0.26 -2.77 -0.15 -2.18 
CO6A2_HUMAN COL6A2 Collagen alpha-2(VI) chain 0.07 0.26 -3.81 -1.96 -3.02 -1.52 
CO6A3_HUMAN COL6A3 Collagen alpha-3(VI) chain 0.12 0.17 -3.11 -2.56 -2.46 -2.01 
CO7A1_HUMAN COL7A1 Collagen alpha-1(VII) chain 0.34 0.10 -1.56 -3.39 -1.20 -2.68 
COPZ2_ HUMAN COPZ2 Coatomer subunit zeta-2 0.06 0.02 -4.07 -5.47 -3.23 -4.37 
CPPED_HUMAN CPPED1 Calcineurin-like phosphoesterase domain- 0.10 0.11 -3.33 -3.15 -2.63 -2.49 
containing protein 1 
D6RJ89_HUMAN ACOX3 Peroxisomal acyl-coenzyme A oxidase 3 0.18 0.12 -2.44 -3.07 -1.92 -2.43 
DCBD2_HUMAN DCBLD2 Discoidin, CUB and LCCL domain- 0.29 0.16 -1.77 -2.68 -1.38 -2.11 
containing protein 2 
DOP2_HUMAN DOPEY2 Protein dopey-2 0.06 0.28 -4.03 -1.82 -3.20 -1.42 
DPYL3_HUMAN DPYSL3 Dihydropyrimidinase-related protein 3 0.10 0.40 -3.27 -1.31 -2.59 -1.00 
E7EUDO0_HUMAN DKK3 Dickkopf-related protein 3 0.17 0.20 -2.57 -2.34 -2.02 -1.83 
ES8L2_HUMAN EPS8L2 Epidermal growth factor receptor kinase 0.34 0.09 -1.56 -3.47 -1.20 -2.75 
substrate 8-like protein 2 
F5GY03_HUMAN SPARC SPARC 0.12 0.15 -3.01 -2.76 -2.38 -2.17 
F8WCU2_HUMAN FKBP7 Peptidyl-prolyl cis-trans isomerase 0.16 0.19 -2.67 -2.38 -2.10 -1.87 
FA49A_HUMAN FAM49A Protein FAM49A 0.13 0.14 -2.96 -2.83 -2.34 -2.23 
FHL1_HUMAN FHL1 Four and a half LIM domains protein 1 0.13 0.09 -2.89 -3.55 -2.28 -2.81 
FHL2_HUMAN FHL2 Four and a half LIM domains protein 2 0.14 0.17 -2.88 -2.58 -2.27 -2.03 
FINC_HUMAN FN1 Fibronectin 0.11 0.13 -3.25 -2.96 -2.57 -2.33 
FKB10_HUMAN FKBP10 Peptidyl-prolyl cis-trans isomerase 0.67 0.07 -0.57 -3.80 -0.40 -3.01 
FKBP10 
FLNC_HUMAN FLNC Filamin-C 0.17 2.65 -2.59 1.41 -2.03 1.19 
FPRP_HUMAN PTGFRN Prostaglandin F2 receptor negative 0.25 0.06 -1.99 -4.09 -1.55 -3.25 
regulator 
(Continued) 
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Table 3. (Continued) 


Swiss-prot 


FUCO_HUMAN 
G3V2M6_HUMAN 


GABT_HUMAN 
GBP1_HUMAN 
GBP2_HUMAN 


GELS_HUMAN 
GLSK_HUMAN 
GPX1_HUMAN 
HOY8D1_HUMAN 
HOYGX7_HUMAN 
H7C2T5_HUMAN 


H7C5L1_HUMAN 
HM13_HUMAN 
HYEP_HUMAN 
ICAM1_HUMAN 
ITA3_HUMAN 
ITAV_HUMAN 
ITPR3_HUMAN 
JAM1_HUMAN 
K1C10_HUMAN 
K1C9_HUMAN 
K2C1_HUMAN 
K7ENN8_HUMAN 
L1CAM_HUMAN 
LAMB1_HUMAN 
LASP1_HUMAN 
LEG3_HUMAN 
LRP1_HUMAN 


MAOX_HUMAN 


MAP1B_HUMAN 
MICA2_HUMAN 


ML12A_HUMAN 
MMP14_HUMAN 
MOT4_HUMAN 
MYH9_HUMAN 
MYL9_HUMAN 
NIBAN_HUMAN 
NMES1_HUMAN 


NXP20_HUMAN 


Gene 
name 
FUCA1 
STAT2 


ABAT 


GBP1 


GBP2 


GSN 

GLS 
GPX1 
PRSS1 
ARHGDIB 
POFUT2 


PTGES2 
HM13 
EPHX1 
ICAM1 
ITGAS 
ITGAV 
ITPR3 
F11R 
KRT10 
KRT9 
KRT1 
TRIM16 
L1CAM 
LAMB1 
LASP1 
LGALS3 
LRP1 


ME1 


MAP1B 
MICAL2 


MYL12A 
MMP14 
SLC16A3 
MYH9 
MYL9 
FAM129A 
NMES1 


FAM114A1 


Protein name 


Tissue alpha-L-fucosidase 


Signal transducer and activator of 
transcription 2 


4-aminobutyrate aminotransferase, 
mitochondrial 


Interferon-induced guanylate-binding 
protein 1 


Interferon-induced guanylate-binding 
protein 2 


Gelsolin 


Glutaminase kidney isoform, mitochondrial 


Glutathione peroxidase 1 
Trypsin-1 
Rho GDP-dissociation inhibitor 2 


GDP-fucose protein O-fucosyltransferase 


2 

Prostaglandin E synthase 2 

Minor histocompatibility antigen H13 
Epoxide hydrolase 1 

Intercellular adhesion molecule 1 
Integrin alpha-3 

Integrin alpha-V 


Inositol 1,4,5-trisphosphate receptor type 3 


Junctional adhesion molecule A 
Keratin, type | cytoskeletal 10 
Keratin, type | cytoskeletal 9 
Keratin, type II cytoskeletal 1 
Tripartite motif-containing protein 16 
Neural cell adhesion molecule L1 
Laminin subunit beta-1 

LIM and SH3 domain protein 1 
Galectin-3 


Prolow-density lipoprotein receptor-related 


protein 1 


NADP-dependent malic enzyme;Malic 
enzyme 


Microtubule-associated protein 1B 


Protein-methionine sulfoxide oxidase 
MICAL2 


Myosin regulatory light chain 12A 
Matrix metalloproteinase-14 
Monocarboxylate transporter 4 
Myosin-9 

Myosin regulatory light polypeptide 9 
Protein Niban 


Normal mucosa of esophagus-specific 
gene 1 protein 


Protein NOXP20 
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M/L 
average 
0.14 
0.19 


0.03 


0.12 


0.09 


0.11 
0.17 
0.74 
0.05 
0.09 
0.13 


1.46 
0.17 
0.26 
0.58 
0.13 
0.12 
0.30 
0.06 
0.13 
0.13 
0.12 
0.21 
0.43 
0.16 
0.16 
0.71 
0.20 


0.31 


0.09 
0.18 


0.09 
0.10 
0.72 
0.13 
0.11 
0.13 
0.40 


0.23 


H/L 
average 
0.16 
0.17 


0.15 


0.13 


0.12 


0.25 
0.37 
0.11 
0.06 
1.91 
0.31 


0.12 
0.26 
0.13 
0.13 
0.19 
0.13 
0.11 
0.08 
0.12 
0.09 
0.11 
0.11 
0.10 
1.20 
0.18 
0.09 
0.12 


0.11 


0.24 
0.26 


0.19 
0.27 
0.06 
0.22 
0.16 
1.19 
0.06 


0.15 


Quantitative Proteome Analysis of Bladder Cells 


log2 M/L 
average 
-2.82 
-2.41 


-5.07 
-3.09 
-3.46 


-3.14 
-2.53 
-0.44 
-4.22 
-3.40 
-2.91 


0.55 
-2.52 
-1.93 
-0.80 
-2.90 
-3.05 
-1.76 
-4.09 
-2.95 
-2.89 
-3.07 
-2.22 
-1.22 
-2.64 
-2.69 
-0.49 
-2.33 


-1.67 


-3.44 
-2.51 


-3.43 
-3.34 
-0.47 
-2.97 
-3.19 
-2.94 
-1.33 


-2.14 


log2 H/L 
average 
-2.60 
-2.52 


-2.72 
72199 
-3.09 


21:99 
-1.43 
-3.14 
-4.01 

0.93 
-1.67 


-3.11 
-1.95 
-2.96 
-2.91 
-2.39 
-2.90 
-3.15 
-3.63 
-3.06 
-3.49 
-3.12 
-3.21 
-3.33 

0.27 
-2.44 
-3.55 
-3.01 


-3.23 


-2.09 
-1.95 


-2.41 
-1.88 
-4.00 
-2.17 
-2.67 

0.25 
-3.95 


-2.73 


z-scores 
M/L 
-2.22 
-1.89 


-4.04 


-2.44 


-2.74 


-2.48 
-1.99 
-0.30 
-3.36 
-2.69 
-2.30 


0.50 
-1.98 
-1.50 
-0.58 
-2.29 
-2.41 
-1.36 
-3.25 
-2.33 
-2.28 
-2.42 
-1.74 
-0.92 
-2.08 
-2.12 
-0.34 
-1.83 


-1.29 


-2.72 
-1.97 


-2.71 
-2.64 
-0.32 
-2.34 
-2.52 
-2.32 
-1.02 


-1.68 


z-scores 
H/L 
-2.05 
-1.98 


-2.14 
-2.31 
-2.44 


-1.55 
-1.10 
-2.48 
-3.18 

0.81 
-1.29 


-2.46 
-1.52 
-2.34 
-2.30 
-1.87 
-2.29 
-2.49 
-2.87 
-2.41 
-2.76 
-2.47 
-2.54 
-2.63 

0.27 
9i 
-2.81 
-2.38 


-2.55 


-1.63 
-1.52 


-1.89 
-1.46 
-3.17 
-1.70 
-2.10 

0.26 
-3.14 


-2.15 
(Continued) 
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Table 3. (Continued) 


Swiss-prot 


OAS3_HUMAN 
PAI2_HUMAN 
PDLI1_HUMAN 
PEA15_HUMAN 
PLSL_HUMAN 
PML_HUMAN 
PNKD_HUMAN 
PPGB_HUMAN 
PTRF_HUMAN 
Q5RHS7_HUMAN 
RAI14_ HUMAN 
RCN1_HUMAN 
RGPS2_HUMAN 


$10A6_HUMAN 
SAMH1_HUMAN 


SELM_HUMAN 
SERPH_HUMAN 
SH3B4_HUMAN 
SNX18_HUMAN 
SNX3_HUMAN 
STC2_HUMAN 
SYCM_HUMAN 


TGM2_HUMAN 


TPM1_HUMAN 
UAP1L_HUMAN 


UBA6_HUMAN 
UBA7_HUMAN 
UN13D_HUMAN 
VAT1_HUMAN 


VIME_HUMAN 
WDFY1_HUMAN 


WIPI1_HUMAN 


Gene 
name 


OAS3 


SERPINB2 


PDLIM1 
PEA15 
LCP1 
PML 
PNKD 
CTSA 
PTRF 
S100A2 
RAI14 
RCN1 
RALGPS2 


S100A6 
SAMHD1 


SELM 
SERPINH1 
SH3BP4 
SNX18 
SNX3 
STC2 
CARS2 


TGM2 


TPM1 
UAP1L1 


UBA6 
UBA7 
UNC13D 
VAT1 


VIM 
WDFY1 


WIPI1 


Protein name 


2-5-oligoadenylate synthase 3 
Plasminogen activator inhibitor 2 
PDZ and LIM domain protein 1 
Astrocytic phosphoprotein PEA-15 


Plastin-2 
Protein PML 


Probable hydrolase PNKD 
Lysosomal protective protein 
Polymerase | and transcript release factor 


Protein S100-A2 
Ankycorbin 
Reticulocalbin-1 


Ras-specific guanine nucleotide-releasing 


factor RalGPS2 
Protein S100-A6 


SAM domain and HD domain-containing 


protein 1 
Selenoprotein M 
Serpin H1 


SH3 domain-binding protein 4 


Sorting nexin-18 
Sorting nexin-3 
Stanniocalcin-2 


Probable cysteine—tRNA ligase, 


mitochondrial 


Protein-glutamine gamma- 
glutamyltransferase 2 


Tropomyosin alpha-1 chain 


UDP-N-acetylhexosamine 
pyrophosphorylase-like protein 1 


Ubiquitin-like modifier-activating enzyme 6 
Ubiquitin-like modifier-activating enzyme 7 
Protein unc-13 homolog D 

Synaptic vesicle membrane protein VAT-1 


homolog 
Vimentin 


WD repeat and FYVE domain-containing 


protein 1 


WD repeat domain phosphoinositide- 
interacting protein 1 


M/L 
average 
0.25 
0.08 
0.15 
0.15 
0.11 
0.14 
0.18 
0.17 
0.15 
1.27 
0.14 
0.46 
0.10 


0.31 
0.35 


0.72 
0.44 
0.11 
0.05 
0.09 
0.24 
0.25 


0.32 


0.06 
0.15 


0.17 
0.15 
0.39 
0.16 


0.15 
0.08 


0.18 


H/L 
average 
0.16 
0.42 
0.39 
0.37 
0.06 
0.15 
0.14 
0.28 
0.53 
0.17 
0.43 
0.16 
0.40 


0.07 
0.11 


0.14 
0.17 
0.08 
1.04 
0.27 
0.17 
0.14 


0.07 


0.15 
0.19 


0.22 
0.08 
0.10 
0.21 


0.38 
0.30 


0.13 


*Proteins shown have at least one z-score value (M/L or H/L) >1.9600 in two biological replicates. 


doi:10.1371/journal.pone.0134727.t003 
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log2 M/L 
average 
-2.02 
-3.64 
-2.76 
-2.74 
-3.20 
-2.79 
-2.48 
-2.58 
-2.74 
0.35 
-2.87 
-1.11 
-3.32 


-1.71 
-1.53 


-0.48 
-1.18 
-3.15 
-4.41 
-3.48 
-2.04 
-1.98 


-1.66 


-4.05 
-2.76 


-2.58 
-2.72 
-1.38 
-2.63 


-2.75 
-3.60 


-2.50 


log2 H/L 
average 
-2.63 
-1.25 
-1.35 
-1.43 
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Fig 4. Functional classification of identified proteins using SWISS-PROT database based on universal GO annotation terms. Proteins shown were 
linked to at least one annotation term within the GO molecular function (A), biological process (B), and cellular component (C) categories. 


doi:10.1371/journal.pone.0134727.g004 
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Confirmation of SILAC and western blotting results by cell staining with 
antibodies 


Antibodies recognizing MAGEA4, THY1, IGF2BP1, VIM, CTNNB1, and FN1 proteins, whose 
expression differed significantly in BC cells vs. HCV29 cells, were used to confirm the results of 
previous analyses and to assess protein distributions. Fluorescence signal intensities in the two 
BC cell lines were significantly higher for MAGEA4, THY1, and IGF2BP1 and significantly 
lower for VIM, CTNNB1, and FN1, in agreement with SILAC, western blotting, and RT-PCR 
results. Preferential localization was observed for MAGEA4 in central cytoplasm (including 
mitochondria and centrosome) and the nuclear region, for THY1 and VIM in cytoplasmic 
membrane and cytoplasm, for IGF2BP1 in the nuclear region and cytoplasm, and for CTNNB1 
and FN1 in cytoplasmic membrane (Fig 7). 


Discussion 


Urothelial carcinoma of the bladder is unique among epithelial carcinomas in its divergent 
pathways of tumorigenesis. At the time of diagnosis for transitional cell BC, ~80% of patients 
are in the low-grade (grade 1-2) non-muscle invasive (cTa or cT1) stage. In patients with low 
grade Ta disease, the 15-year progression-free survival is 95% with no cancer-specific mortality 
[36]. Therefore, biomarkers are needed for predicting a risk of stage progression from non- 
invasive to invasive, or non-metastatic to metastatic and predicting responsiveness to systemic 
therapies. The human cell lines HCV29 (normal bladder epithelia), KK47 (low grade nonmus- 
cle invasive bladder cancer), and YTS1 (metastatic bladder cancer) have been widely used in 
studies of molecular mechanisms and cell signaling during the progression of bladder cancer to 
muscle or metastatic states [37]. However, little attention has been paid to global quantitative 
proteome analysis of these three cell lines. 

Quantitative analysis of proteins at different stages of BC progression is a challenging but 
important task for understanding of disease mechanisms. SILAC, a differential isotope labeling 
strategy that involves metabolic labeling of proteins in vivo, has been widely applied in cell biol- 
ogy and studies of model organisms such as yeast, bacteria, nematodes, plants, and mice [38, 
39]. In SILAC, natural “light” isotopes of carbon, nitrogen, and hydrogen in amino acids incor- 
porated during protein translation are substituted with "heavy" isotopes such as 13, 15N, and 
°H. No study to date has quantified differences in protein abundance at various stages of BC. 
Previous studies have focused on comparative protein levels in BC patients vs. control subjects, 
but not on alterations of protein levels during BC progression [8,9]. 

Identification and characterization of protein levels at various steps of differentiation are 
essential for our understanding of normal tissue development and malignant transformation. 
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Fig 5. Functional network analysis of differentially regulated proteins with z-score cutoff 95% observed in different stages of BC cells using 
Ingenuity Pathways Analysis (IPA). (A and B) Top network functions of DNA replication, molecular transport, cell growth, and cell proliferation for 
upregulated proteins. (C and D) Top network functions of cellular movement, immune cell trafficking, and lipid metabolism. Solid lines: direct known 
interactions. Dashed lines: suspected or indirect interactions. White: proteins known to be in the network but not identified in our study. 


doi:10.1371/journal.pone.0134727.g005 


In the present study, we applied the SILAC method to analysis of HCV29, KK47, and YTS1. 
This strategy provided protein information for all three cell lines during one MS experiment. 
Following normalization by the z-score method, differential regulation was observed for 110 
proteins in KK47 vs. HCV29 and for 87 proteins in YTS1 vs. HCV29. These differentially regu- 
lated proteins, which may play important roles in BC development, include SFN (14-3-3 pro- 
tein sigma), SELENBP1 (selenium-binding protein 1), COL6A3 (collagen 03 (VI) chain), and 
CD70. COL6A2 and COL6A3 protein levels are reduced in urine of BC patients [40]. SFN is 
downregulated in invasive bladder transitional cell carcinomas undergoing epithelial-to-mes- 
enchymal conversion and highly upregulated in pure squamous cell carcinomas [41]. Some of 
the proteins that are altered in BC are also altered in other cancers. In particular, expression of 
FOLRI (folate receptor alpha) is increased relative to normal tissue in many types of epithelial 
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Fig 6. Confirmation of SILAC-determined protein and RNA abundances. (A) SILAC L:M:H average ratios for six selected proteins. (B) Western blotting 
analysis of selected proteins. Proteins were transferred to a PVDF membrane, probed with their primary antibodies, and incubated with HRP-conjugated 
rabbit anti-mouse or goat anti-rabbit secondary antibodies. (C) Densitometric quantitation of the protein levels. The protein was normalized by the tublin, and 
then compared to HCV29, which were arbitrarily set at 1.0. (D) Gene expression for the proteins was analyzed by qRT-PCR. Relative expression in 
comparison to control samples was analyzed by the 2-A4Ct method and represented as Logs. Expression of genes above Logo(2) or below Logo(1/2) was 
significantly upregulated or downregulated, respectively. 


doi:10.1371/journal.pone.0134727.g006 


cancer, including non-mucinous ovarian, endometrial, non-small cell lung, colorectal, and 
breast cancers [42]. 

The DAVID bioinformatic resources provide a set of powerful tools to explore their large 
gene lists in depth from many different biological perspectives in order to fully extract associ- 
ated biological meanings. Enrichment analysis of GO terms in our identified proteins indicated 
that neutral amino acid transmembrane transporter activity and metabolic processes were sig- 
nificantly upregulated and that protein binding, cell adhesion, biological adhesion, and 
immune response were significantly downregulated in the BC cell lines. The observed changes 
in GO terms tended to promote cellular proliferation, tumor development, cancer cell progres- 
sion and metastasis, and escape from immune system surveillance. 

A major challenge in cancer biology is the formulation of biological hypotheses regarding 
biomarker candidates [43]. Pathway analysis of differentially expressed proteins indicates that 
DNA replication and molecular transport, cell growth, and cell proliferation are requirements 
for cancer cell metastasis. Intercellular adhesion molecule-1 (iCAM-1) is a transmembrane gly- 
coprotein present at basal levels in a wide variety of cell types and is upregulated in response to 
a number of inflammatory mediators [44]. The biological significance of iCAM-1 expression in 
cancer remains controversial; it is elevated in gastric, breast, oral, and thyroid cancer tissues 
[44-47] but reduced in some ovarian adenocarcinoma cell lines and primary tumors [48]. 
Also, iCAM-1 expression is up-regulated in squamous cell types associated with inflammation, 
such as schistosomal bladder cancer [49]. However, in the present study, iCAM-1 expression 
was reduced in metastatic bladder cancer cells. Moreover, CIT, a novel tissue-specific Ser/Thr 
kinase that encompasses the Rho-Rac-binding protein Citron, plays a role in cytokinesis and in 
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Fig 7. Differential expression revealed by cell staining. HCV29, KK47, and YTS1 cells were cultured and 
stained with six antibodies directed to identified proteins (MAGEA4, THY1, IGF2BP1, VIM, CTNNB1, FN1) 
labeled with Cy3 as described in M&M. Images are shown of merge images of Cy3-conjugated antibodies 
and DAPI staining of nuclei (objective magnification 60x). Scale bars: 70 um. 


doi:10.1371/journal.pone.0134727.g007 
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Rho signaling that modulates myosin phosphorylation and cell adhesion [50]. Upregulation of 
CIT promotes DNA replication and cancer cell proliferation. 

In conclusion, we successfully applied the SILAC method to identify and quantify proteins 
whose level is significantly up- or downregulated during BC development or progression. Simi- 
lar advanced proteome techniques will be useful for further elucidation of biomarkers and 
molecular mechanisms in BC and other types of cancer. 
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